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ABSTRACT 

Structural models are examined for the influence of a ring with an asymmetrical 
cross section on the linear elastic response of an orthogonally stiffened cylindrical 
shell subjected to internal pressure. The first structural model employs classical 
theory for the shell and stiffeners. The second model employs transverse shear 
deformation theories for the shell and stringer, and classical theory for the ring. 
Closed-end pressure vessel effects are included. Interacting fine load intensities 
are computed in the stiffener-to-skin joints for an example problem having the 
dimensions of the fuselage of a large transport aircraft. Classical structural theory 
is found to exaggerate the asymmetric response compared to the transverse shear 
deformation theory. 

INTRODUCTION 

The cabin pressurization in a transport aircraft causes about a 10 psi pressure 
differential across the skin. An unstiffened, or a monocoque fuselage would carry 
this internal pressure load as a shell in membrane response, like a pressure ves- 
sel. However, internal longitudinal and transverse stiffeners are necessary to carry 
loads due to flight maneuvers, landing and ground handling, etc. How the loads 
are transferred in the stiffener-to-skin joints under pressurization is necessary for 
determining the load capacity of these joints. 

* Graduate Research Assistant, Aerospace and Ocean Engineering 
f Professor of Aerospace and Ocean Engineering 
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The design of stiffener-to-skin joints was cited by Jackson, et al. (1984) as 
one of the major technology issues in utilizing graphite/epoxy composites in the 
fuselage of a large transport aircraft. Stiffeners can be attached to the skin by 
either fasteners, co-curing, adhesive bonding, or some combination of these methods. 
Where fasteners are required in a graphite/epoxy structure, aluminium fasteners 
cannot be used because of galvanic corrosion to the metal. More expensive fasteners, 
like titanium, are required to avoid corrosion. 

A ring, or frame, with asymmetrical cross section is commonly used as a trans- 
verse stiffening member in the fuselage of a transport aircraft. The influence of this 
ring asymmetry on the distribution of the interacting loads in the stiffener-to-skin 
joint is the subject of this paper. Two structural models are considered. The first 
model employs classical theory for the shell and the stiffeners. The second model 
employs transverse shear deformation theories for the shell and stringer, and classi- 
cal theory for the ring. We have previously published results for the response with 
symmetrical section rings and stringers using classical theory (Johnson and Rastogi, 
1994). 

Hence, the objective of this paper is to examine structural models for the linear 
elastic response of an orthogonally stiffened, composite material cylindrical shell 
subjected to internal pressure, where the ring has an asymmetrical cross section 
and the stringer has a symmetrical cross section. This objective is part of a larger 
effort to develop an analysis/design capability for the stiffener-to-skin joint of a 
large transport aircraft. A potential benefit of such an analysis/design capability is 
to use fewer expensive fasteners in the graphite/epoxy fuselage. 

MATHEMATICAL MODEL 

An idealized mathematical model is assumed for the semi-monocoque fuselage 
to study the generic characteristics of the response in the vicinity of the stiffeners’ 
intersection. The model is of a very long circular cylindrical shell internally stiffened 
by identical stringers equally spaced around the circumference, and identical frames 
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or rings, equally spaced along the length. In general, the spacing of the stringers is 
not the same as that of the rings. The structure is periodic both longitudinally and 
circumferentially, and the loading is spatially uniform. Consequently, a structural 
repeating unit can be defined whose deformation determines the deformation of the 
entire structure. A typical repeating unit consists of a portion of the shell wall 
centered over the portions of stringer and ring as shown in Fig. 1. The radius of 
the middle surface of the undeformed cylindrical shell is denoted by R , and the 
thickness of the shell is denoted by t. Axial coordinate x and the circumferential 
an gle 0 are lines of curvature on the middle surface, and the thickness coordinate 
is denoted by z, with —tj 2 < 2 < t/2. The origin of the surface coordinates is 
centered over the stiffeners intersection so that —l<x<l and —0 < 6 < 0, where 
21 is the axial length, and 2RQ is the circumferential axe length of the repeating 
unit. 

The stiffeners are mathematically modeled as one-dimensional elements, or dis- 
crete beams, so that the actions transmitted by the stiffeners to the inside of the 
shell wall are represented by distributed line load intensities. In this paper it is 
assumed that the stringer is symmetric about the x-z plane through its centroidal 
axis and the ring is asymmetric. On the basis of the symmetry about the x-axis 
for the unit, only the interacting line load components tangent and normal to the 
stringer are included in the analysis. However, due to an asymmetric ring, the 
interacting line loads between shell and the ring consist of three distributed force 
components and a tangential distributed moment component. The shell-stringer 
interacting force components per unit length along the contact lines are denoted 
by A IS (x) for the component tangent to the stringer and A 2S (x) for the component 
normal to the stringer. The three shell-ring interacting force components per unit 
length along the contact lines axe denoted by A ir (0) for the component acting in 
the axial direction, A# r (0) for the component tangent to the ring, and \ zr (0) for 
the component normal to the ring. The shell-ring interacting moment component, 
tangent to the ring, per unit length along the contact line is denoted by Ae r (0). 
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These interacting loads acting in a positive sense on the inside surface of the shell 
are shown in Fig. 2. The purpose of the analysis is to determine these distributed 
line load intensities and also, to examine the differences in their distribution for the 
two structural models described earlier. 

For both models, the linear elastic response of the repeating unit to internal 
pressure is obtained by utilizing Ritz method and the principle of virtual work 
applied separately to the shell, stringer, and ring. The virtual work functionals are 
augmented by Lagrange multipliers to enforce kinematic constraints between the 
structural components of the repeating unit. The Lagrange multipliers represent 
the interacting line loads between the stiffeners and the shell. Displacements are 
separately assumed for the shell, stringer, and the ring. 

TRANSVERSE SHEAR DEFORMATION FORMULATIONS 
Shell 

A consistent first order transverse shear deformation theory is developed to 
model the shell. Based on the assumption that the shell thickness t is relatively 
small and hence, does not change during loading, the displacements at an arbitrary 
material point in the shell are approximated by 


U(x, 0, z) = u(x, 6) + z<f> x (x,0) (1) 

V(x, 0, z) = v(x, 9) + z<f>s(x, 6) (2) 

W(x,e,z) = w(x,0) (3) 


where u(x,0), v(x, 0) and w(x,6) are the displacements of the points of the ref- 
erence surface, and 0 £ (x, 0) and Oo( x . 9) are the rotations of the normal to the 
reference surface as shown in Fig. 3(a). Using Eqs. (1) to (3) and assuming small 
displacement gradients, the three-dimensional engineering strains are 


— ^7 


+ 


C00 


£06 + ZKqq 

(1 + *) 


= 0 (4) 
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The transverse shear strains e xz and eg z represent average transverse shearing 
strains through the thickness of the shell since Eqs. (1) to (3) are approximate 
in the z-coordinate. In Eqs. (4) to (6), the two-dimensional, or shell, strain mea- 
sures, which are independent of the z-coordinate, are defined by 


du 

_ d(j) x 


^xx ~ Ti 
OX 

^XX — 

ox 

(0 

1 dv w 

1 d(f>e 

(8) 

e " = Rd0 + R 

K " ~ R dO 


dv 1 du 

(9) 

7x9 = 

dx + RdO 

H 

II 

Cb 

N 

1^ 

1 d(f> x 1 dv 

+ R dO + Rdx 

(10) 

d(f>e 

^xO — o 
OX 

1 d<f> x 1 dv 

R dO Rdx 

(11) 

dw 

, v 1 dw 

(12) 

~ixz == 4*x H 

70z = H ~ r + 


If we set the (average) transverse shear strains in Eq. (6) to zero, then the 


rotations of the normal are 


dw 

Yx — o 

OX 

v 1 dw 
^ e=Z R~ R~dO 


(13) 

(14) 


so that 


K'xO — K'xO — 


2 d 2 w 2 dv 
R dxdO R dx 


k x g = 0 


(15) 


Hence, the thickness distribution of the shear strain reduces to 


^xO 


7x0 + 2(1 + 2 

(! + i) 


( 16 ) 


which coincides with the results of Novozhilov’s (1964) classical shell theory. 



It is evident from Eq. (5) that three shell strain measures are needed to rep- 
resent the shear strain distribution through the thickness in the transverse shear 
deformation shell theory. Whereas, only two shell strain measures are required in 
classical shell theory to represent the shearing strain distribution through the thick- 
ness (refer to Eq. (16)). Also it can be shown that under rigid body rotation of the 
shell, the nine shell strain measures, given by Eqs. (7) through (12) vanish. (For 
Novozhilov’s classical shell theory, six shell strain measures given by Eqs. (7-9) and 
(15) vanish under rigid body rotations). 

The physical shell stress resultants and stress couples in terms of stress compo- 
nents are given, in usual way, by 


(N xx , M X x) — J (1) z)(T xx (l + j^) dz 
(Nee, Mee) = J(l, z)a 0 e dz 
(N x0 , M x0 ) =j(l,z)a x e(l + ^)dz 
(Ne x , M Bx ) =£ (l,z)a 0x dz 

Qx = ^ & X z( 1 *b ^) dz 

Qe =Ja Sz dz 

A generalized 9x1 stress vector for the shell is defined by 


(17) 


& shell — [AT**, Neg, N 0X ,M XX ,M 00 , M x $, M x0 , Qx, Qe] 


(18) 


in which M x0 and M x e are the mathematical quantities conjugate to the modified 
twisting measures R x e and k x0 , respectively, and are defined in terms of the physical 
stress couples by 

M x0 = \(M xe + Me x ) M x0 = \(M xS - M 0X ) (19) 
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The nine elements of the stress vector in Eq. (18) and the relations of Eq. (19) 
determine all the stress resultants and stress couples listed in Eq. (17) except for 
shear resultant N x q. The shear stress resultant N x $ is determined from moment 
equilibrium about the normal for an element of the shell. This so-called sixth 
equilibrium equation is 


n z6 = n 6x + 


M $x 

R 


( 20 ) 


The generalized strain vector for the shell is 

^ . rp 

£ shell = 


( 21 ) 


This strain vector is conjugate to the stress vector in the sense that the internal 
virtual work for the shell is given by 


shell ~ j j s ^ ‘shell ® shell dS 


( 22 ) 


where S denotes the area of the reference surface and dS = dxRdd. This expression 
for the internal virtual work can be derived from three-dimensional elasticity theory 
by using Eqs. (4) to (6) for the thickness distributions of the strains and the 
definitions of the resultants given by Eqs. (17) and (19). 

Consistent with the transverse shear deformation theory, the linear elastic con- 
stitutive law for a laminated composite shell wall is given by 
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and 


f Qx 1 _ A 44 ^45 1 f 7 12 1 

\ Qo J . -d-45 ^55 J \ 7 6z J 


(24) 
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in which stiffnesses Ai } ,B l} and D tJ are given in Appendix. The transverse shear 
stiffnesses, A 44 , A 45 , and A 55 can be calculated by two different methods. The first 
method is based on the assumption of constant transverse shear strain distribu- 
tion through the thickness, and the second method is based on the assumption of 
constant transverse shear stress distribution through the thickness. In the present 
analysis, we have used the first method to compute the transverse shear stiffnesses. 

The statement of virtual work is 


sw 


int 

shell 


= sw; xt + SW{ xi 


(25) 


where the external virtual work for a cylindrical shell under constant internal pres- 
sure, including an axial load due to the closed-end effect, is written as 


© 

SWp Xt = JJ p Sw dS + p J ^-dO [8u(l,0) - 8u(-l,0)] 
' 5 -© 


(26) 


and the external (or augmented) virtual work due to the interacting loads is 
/ 

SWl xt = j |A XJ (x)[6u(x 1 0) - -6cf> x (x,0)\ + A^ a (x)6u;(x,0)| dx 

-i 

© 

+ J {\ X r(e)[Su(0 : 9)- t -SM^0)]+\er(0)[Sv(0,9)- t -SM^0)^ 

-0 

j 

+ \ 2T (0)6w(0,e) - Ag r (0)f>( ^ | r=(J )} (R- ^)<W - Q[hu(l,(>) - «u(-/.0)] 

(27) 

The axial force Q in Eq. (27) is an additional Lagrange multiplier that accounts 
for axial load sharing between the stringer and shell. 


Stringer 

Stringer displacements u 3 (x ) and w s (x), and the rotation of the normal <f>e 9 (x) 
are shown in Fig. 3(b). Based on transverse shear deformation theory, the virtual 
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work expression for the stringer is 


l 

j "h “f h 5 8 *j zs \dx 

-l 


l 

J {w x)[Su 3 (x) + e s 8 <t>e s (x)\ 

-l 


+ A zs (x)^tx; a (x)| dx + Q[ 6 u s (l ) - Su s (-l)] 

(28) 

in which N xs is the axial force in the stringer, Mg s is the bending moment, V a is the 
transverse shear force, e xs is the normal strain of the centroidal line, the product 
ZK 0 S is the portion of the axial normal strain due to bending, 7 zs is the transverse 
shear strain, and e s is the radial distance from the stringer centroid to the contact 
line along the shell inside surface. The strain-displacement relations and Hooke’s 
law for the stringer are 


= u 


K 6s — 4>' 6 s 


' )zs — &9s T W s 


(29) 


N X3 = ( EA) s e xs Mes = ( EI) 3 k 9s V s = (GA) slzs (30) 

in which the prime denotes an ordinary derivative with respect to x. 

CLASSICAL FORMULATIONS 


Shell 

The shell is modeled with Sanders’ (1959) theory for thin shells. Define a gen- 
eralized strain vector in terms of the shell strain measures by 

t 3 hell = [^xxi £ 661 'Yxdi ^111 K' 661 ^x#] (31) 


The first five strain measures of the shell reference surface in Eq. (31) are related 
to the displacements by Eqs. (7-9), and the sixth strain measure, k x o, is given by 

__ d<t>e , 1 d<f> x 1 

+ r W + r* 



where the rotation about the normal, <f>, is 

1 dv 1 du 

* = ~ HHd’ 


( 33 ) 
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and the rotations <j> x and <j>g of the normal are given by Eqs. (13) and (14). 

Define a generalized stress vector in terms of the stress resultants and couples 
of Sanders’ theory by 

v shell — [-^m ft 66, M xx , Mgg, M*g] T (34) 

such that the internal virtual work is given by Eq. (22). Quantities N* e and M s x6 are 
the modified shear resultant and twisting moment resultant in the Sanders theory. 
Hooke’s law for a laminated composite shell wall is 

v shell = He s hell H — g T p (35) 

in which the 3 x 3 sub-matrices A, B and D are given by classical lamination theory 
(Jones, 1975). The external virtual work expressions for the classical shell theory 
are still given by Eqs. (26) and (27), but the rotations in Eqs. (27) are given by 
Eqs. (13) and (14). 

Stringer 

The stringer is modeled with Euler-Bernoulli beam theory thereby neglecting 
the transverse shear strain. Hence, equating y zs in Eq. (29) to zero results in the 
following expression for (f>g s . 

4>9s = -w', ( 36 ) 

It may be noted that neglecting the transverse shear strain would also modify the 
virtual work statement given by Eq. (28), and the third equation in the Hooke’s 
law, Eq. (30), is neglected. 

Ring 

Ring displacements are denoted u r (0), v r (@), and u> r (0), and the rotations are 
denoted by (f> xr (Q), <f>e r (0), and 4> zr (8) as shown in Fig. 3(c). The structural model 
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is based on Euler- Bernoulli hypotheses. The statement of virtual work is 


0 

J g r Seg r T M xr 8 k xt T M zr SK zr + T r 8T r \Ro d6 
-© 


© 


- J | A X r(d)[8u r (8) + e r 8cj)o r (d)] + A $ r (9)[8v r (0) + e r 6<f> xr (9)] 
-© 


( 37 ) 


+ 


A zr (0)Sw r (0) + yW(0)^*r(0)} (1 + jjr) d6 


in which N# r is the circumferential force, M xr is the in-plane bending moment, 
M zr is the out-of-plane bending moment, T r is the torque, tg r is the circumferential 
normal strain of the centroidal arc, n xr is the change in curvature due to in-plane 
bending, K zr is the change in curvature due to out-of-plane bending, r T is the twist 
rate, e r is the distance from the ring reference arc to the contact line along the 
shell inside surface, and R 0 is the radius of ring reference arc. The rotations and 
strain-displacement relations are 


€0r = ^ r + Wr) 

Kq 

T r — ~^~{4 > 0r T <f>zr) 
Kq 


^ir jy 

Kq 

fixr ~ rs (^r d) r ) 
Kq 


tr = - <f>0r) 

Kq 

(j>zr = d r 

Kq 


(38) 


in which (j) xr is the rotation around x-axis, <f> zr is the rotation around z-axis, and 
the over-dot denotes an ordinary derivative with respect to 9. Hooke’s law is 

Ng r = (EA) r €g r M X r = (^n)r^rr {El zx )r K ‘ zr 

M zr = {El zz ) r K zr — ( EI zx )r K xr T r = ( GJ) r T r 


(39) 


DISPLACEMENT CONTINUITY 

In order to maintain continuous deformation between the inside surface of the 
shell and stiffeners along their lines of contact, the following displacement continuity 
constraints are imposed: 

Along the shell - stringer interface (i.e., —l<x<l,0 = 0), 

g xs = u(x, 0) - ^ <f> x (x , 0) - [u 3 (x) + e s <f>e s (x)] = 0 (40) 
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gzs = w(x,0) - w a (x) = 0 (41) 

Along the shell - ring interface (i.e., x = 0, -0 < 0 < 0), 

g xr = u(0, 6) - ^<t>x( 0, 0) - [u r (fl) + e r <t>or(0)] = 0 (42) 

g 9r = u(0, e) - *-M 0. °) - [t>r(0) + e r <f>xr(0)\ =0 (43) 

g ir = w(0,6)-w r {6) = 0 (44) 

«.- = -tL=,-W‘> = 0 (45) 

The variational form of these constraints are 

l 

f [S\ x3 g xs + fi^zs9zs] dx = 0 (46) 


J \6\ xr g X r + fi^OrgBr + fi^zrQzr + SAorGor] (Ro + e r) dd = 0 (47) 

-0 

The constraint that the elongation of the shell at 6 = 0 and the elongation of 
the stringer are the same is 

6Q{[u(l,0)-u(-l,0)]-[u 9 (l)-u s (-l)}} =0 (48) 


DISPLACEMENTS, ROTATIONS, AND INTERACTING 
LOAD APPROXIMATIONS 


The periodic portions of the displacements and rotations are represented by 
truncated Fourier Series having fundamental periods in the stringer and ring spac- 
jjijr The non-periodic portions of the displacements due to axial stretching are 
represented by simple terms in x. The Fourier series reflect symmetry about the 
x-axis for the repeating unit. For the shell, displacements of the middle surface (see 
Fig. 3a) are represented as 


M N 


M N 


m n Sin(a m x)Cos(f3 n 8)-\- y] ^ ^ U 2 mnGos(amX)Cos(f 3 n 8 ) 

(49) 


m— 1 n=0 


m=l n=l 
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M N 


M N 


.(«.*) = EE EE V2mnSin(at m x)Sin(fi n d ) (50) 


m = 0 n=l 

M N 


m=l n=l 

M N 


w(x,0)= ^2 ^2 w imnCos(a m x)Cos(l3 n 0) + ^ y: w 2mn 5m(Q m x)Co^(/9 n 6>) 

(51) 


771=0 77=0 


m=l n=l 


and rotations of the normal are 
M N 


M N 


= EE*>- n 5m(a m x)C'os(/3„^) + ^2 4>x2mnCos(a m x)Cos(f3„0) 

771 = 1 ri—0 m=l 77=1 

(52) 

M N M N 

<f) 6 (x,d)= ^>2 ^2<l>9imnCos(amx)Sin(P n 8)+ ^2 Y2 <t>62m n Sin(a rn x)Sin{0 n $) 

m=0 77 = 1 m — 1 77=1 

(53) 

in which a m = and /?„ = ~ where m and n are non-negative integers. Note 
that some terms in the truncated Fourier Series of Eqs. (49-53) have been omitted. 
The coefficients of the omitted terms are U 200 , u 2m o, «20n, »2m0, 4>x 200 , <f>x2mo, and 
4 > x20ni in which m — 1,2 , M and n — 1,2, ..., N. The rationale for their omission 
is discussed in the following sub-section. The displacements of the centroidal line 
of stringer (see Fig. 3b) are 
M 


M 


uj,x) = u 3lm Sin(a m x) + ^ u 3 2 mCos(a m x ) 


(54) 


777 = 1 


777 = 1 


M M 

w s (x) ~ E IV m StTX^Ot fji X ) ( ^V s2 tn CoS ( Oi X ) (55) 

m = 1 m= 1 

and the rotation of the normal of the stringer about the 0-axis is 

M M 

4>ds{x) = <t> 0 slm Sin(a m x ) + '22 < f > 0 s 2 mCos(a m x ) (56) 

777 = 1 777 = 1 

where the coeffcients u s 20 , w S 2 o and <f> 0 s 2 Q ar e omitted. Coefficient qo in the axial 
displacement field of the shell and q\ in the axial displacement field of the stringer 
represent elongations of each respective element caused by either an axial mechanical 
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load or due to close-end pressure vessel effects. The displacements of the reference 
circle of the ring (see Fig. 3c) are 


N 

M«) = E “ rn Cos(/3 n 9 ) 

71=1 


(57) 


and twist of the ring is 


N 

v r {9) = ^ v rn Sin(0 n 9 ) 

n = 1 

N 

W r (9) = W rn Cos(l3 n 6 ) 

71=0 


N 

<t>er(9) = <f>8rnCos(0 n O) 

n=l 


(58) 

(59) 


(60) 


where the coeffcients u r 0 and <p$ r o are omitted. The distributions of the interacting 
loads, or Lagrange multipliers, are taken as 


M m 

Ajj(l) = y ^ ^xslmSlTl(^Ol m X^ T y ^ ^x32mCos(cXmX) (61) 

m = l m= 1 

M M 

^zs(x) = ^ ^zsl m SlTli^CX xnX') 4“ ^ ^zs2m Cos(a m x) (62) 

m=l m = 1 

N 

M») = E ^xrnCos(f3 n 9) (63) 

n = l 
N 

\ 9r ($) = Y i \ 9 mSin(0 n 6) (64) 

n=l 

N 

M«) = E A zrn Cos{M) (65) 

71=0 

N 

Aor{9) = £ ^rnCostM) (66) 

n=l 

where the coefficients A X 2sO> A X 2 «o> A xr o, and T^ r o are omitted. 
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Terms Omitted in the Fourier Series 

Terms omitted in the truncated Fourier Series for the displacements, rotations, 
and the interacting loads were determined from rigid body equilibrium conditions 
for the ring and stringer, and from displacement continuity conditions between the 
shell and the stiffeners. The external virtual work for the stringer and ring must 
vanish for any possible rigid body motions of these elements. For the stringer these 
rigid body motions are spatially uniform z-direction and z-direction displacements. 
(A rigid body rotation of the stringer in the z-z plane is not considered since this 
motion would violate longitudinal periodicity of the repeating units.) Vanishing of 
the external virtual work for an arbitrary rigid body displacement of the stringer 
in the axial direction leads to the z-direction equilibrium equation 

/ 

J A IS (z ) dx = 0 (67) 

-l 

Similarly, the equilibrium equation for a rigid body displacement of the stringer in 
the z-direction is 

J A zs (z) dx = 0 (68) 

-z 

If the ring is considered in its entirety, that is, as made up of an integer number 
of repeating units around its circumference, the rigid body motions that lead to non- 
trivial equilibrium conditions are a displacement in the z-direction and a rotation 
about the z-axis. The z-direction equilibrium equation is 

© 

J A xr (9) ( Rq + e r ) dS = 0 (69) 

-0 

and the moment equilibrium equation about the z-axis is 

© 

J X er {6) ( R 0 + trf dd = 0 (70) 

-© 
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Equilibrium Eqs. (67) to (69) imply that coefficients 


A I3 2o = 0 A 2S 2o — 0 A rr0 — 0 (71) 


in the Fourier Series for the interacting loads, and these conditions have been taken 
into account in Eqs. (61) to (63). The sine series for \ Sr given in Eq. (64) satisfies 
the equilibrium condition given in Eq. (70). 

Consider the variational form of the constraints, Eqs. (46) and (47), for the 
constant components of the virtual interacting loads. These equations are 


N 


w 200 — z200 T ^ ^ (^20n — ^ < t > x20n ) ( u s20 T e s4’9a2o) 


S A 


xs 20 


= o 


(72) 


n=l 


N 



W\0n 


- Ws20 


8 A 


zs 20 


0 


(73) 


M 


u 200 ~~ 200 + y^( w 2m0 — ~<f ) x2mo) ( w r0 + e r ( l ) 0r0 ) <^xr0 0 

m=l 


(74) 


Since these equations axe satisfied on the basis that ^A I3 20 — 0? &^zs20 — 0 and 
S\ xr0 = 0, consistent with Eq. (71), the bracketed terms in Eqs. (72) to (74) do not 
necessarily vanish. The implication that these bracketed terms in Eqs. (72) to (74) 
do not vanish is that displacement continuity conditions axe not satisfied pointwise. 
Pointwise continuity can be achieved by taking each Fourier coefficient appearing 
in the bracketed terms of Eqs. (72) to (74) to be individually zero. Fourier Series 
given in Eqs. (49), (52), (54), (56), (57), and (60) reflect this choice. Moreover, 
Fourier coefficients u 2 0 0 . U32O1 and u r 0 represent rigid body displacement in the axial 
direction for the shell, stringer, and ring, respectively, and setting them to zero can 
be justified on the basis that rigid body displacement does not contribute to the 
deformation of the structural elements. Since Fourier coefficient w s 20 repersents 
rigid body displacement of the stringer in the 2 -direction, it would seem that it 
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should be set to zero as well. However, to maintain continuity between the stringer 
and the shell in the z-direction, we impose the condition 


N 

^2 W 10n-Ws20 = 0 ( 7 5 ) 

n=0 

to determine w s 20 after obtaining the solution for the displacement components 
that deform the shell; i.e., Fourier coefficients uq 0 ni « = are taken to be 

non-zero independent degrees of freedom since the stringer coefficient w s2 0 is not a 
part of the solution vector. 

Finally, consider the constraint equation associated with 8Ag r o, the spatially 
uniform component of the interacting moment intensity, which was omitted in the 
series given by Eq. (66). Derived from Eq. (47), this constraint equation is 



rn = 1 


w 2m0 + </>0rO 


8Ag r0 — 0 


(76) 


We equated the constant component of the twist, <j>g r0 , to zero from the consider- 
ations associated with Eq. (74). Consequently, a non- zero value of the constant 
component of the interaction moment intensity, Ag r 0 ^ 0, would not contribute 
to the equilibrium of the ring, since Ag r o and <j)g r o are conjugate variables in the 
external work for the ring (refer to Eq. (37)). Since <f>g r 0 = 0, it is necessary that 
Ag r o = 0 to achieve consistent conditions for the torsional and out-of-plane bending 
equilibrium of the ring. With 8Ag r o = 0 in Eq. (76), the bracketed term does not 
necessarily vanish, and as a result pointwise rotational continuity betwen the shell 
and the ring is not assured. Pointwise rotational continuity is achieved if we take 
the coefficients w 2m o = 0, m = 1 , as was done in the Fourier Series for the 

normal displacement of the shell given by Eq. (51). 


DISCRETE EQUATIONS AND THEIR SOLUTION 


Transverse Shear Deformation Model 
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The discrete displacement vector for the shell is the (10MN + 3M + 3A r + 2) x 1 


vector 


in which subvectors are 


T i T 


r ~T - T »Tl 

shell — [«0 5^1 


(77) 


u 0 = [9o,^ioo,uioi,*«ioi,^ioi,--- 5 wioiV,^ioAr,^ioN] 


(78) 


Um — [wimO y ^lmO 5 ^ilmO y^lml y ^ 2 mli ^lmli ^2 ml 5 ^lml 5 ^2 ml 1 ^i2mh 

^lmli ^2ml ? •••? w 2mN y ^xlmN y (79) 

$x2mN 1 $6\mN s $02mN\ 

where m = 1, ..., M 

The (6M + 1) x 1 discrete displacement vector for the stringer and (4 N + 1) x 1 
vector for the ring are 

'll str — [ijl y H si 1 y H 921 y sll y m s21 y $ 0 sll y $0921 y ••• y H si M y H s2 M y 


WslMyWs2My<POslM 


y $0 s2m\ 


u 


ring 


= [ W r o,U r i,V r i,W r i,$0 r i,...,U r NyV r SyWrNy$OrN ] 


(80) 


(81) 


in which the term ie„o for the stringer has been omitted as discussed in reference 
to Eq. (75). The 4 M x 1 discrete interacting loads vector for the shell-stringer 
interface and (AN + 1) x 1 vector for the shell-ring interface axe 

A sir = si 1 y ^xs21 y ^ zsl 1 y ^ zs2\ y • • « y si M y ^xs2 M y M y ^zs2 M ] (^2) 


^rtn^ — zrO y ^xrl y y ^ zrl y ^Orl y • • • y ^xrN y ^OrN y ^ zrN y -^Or 2v] 


(83) 


Classical Model 

The discrete displacement vector for the shell is the ( 6MN + 2 M + 2N + 2) x 1 
vector 

V.T „,T ,T1^ (84) 


r ~T - T -T 1 

9kell [«0 9^1 9 ■•*9 ^A# J 
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in which subvectors are 


u o 


= [90 , u>ioo, uioi,u>ioi,...,wion,wion] 


(85) 


U T 


= [uimO 1 w lmO ) u lm\i u 2mli y lml) ^ml j w lml 1 w 2ml 1 •••) u lraN r u 2mN 1 


^lraJVi^mNi^lmJV 


, W2mN 


( 86 ) 


where m = 1 , M 

The (4M + 1) x 1 discrete displacement vector for the stringer and (477 + 1) x 1 
vector for the ring are 


Ustr = [ qi,U s n,U s 21,W 3 n,W s 2\,...,U s iM,U s 2M,U)slM,W s 2M ] (87) 


taring 


= \w r0 ,u rl , V r i , U+l , <j>Br 1, •••, «rN,«rJV, l"rN, <?Wiv] 


( 88 ) 


The 4M x 1 discrete interacting loads vector for the shell-stringer interface and (47V+ 
1) x 1 vector for the shell-ring interface are the same as for the shear deformation 
model and are given by Eqs. (82) and (83). 

The approximations in Eqs. (49) through (60) for the displacements and Eqs. 
(61) through (66) for the interacting loads are substituted into the virtual work 
functionals for each structural element, and also substituted into the variational 
form of displacement continuity constraints. Then integration over the space is 
performed. (The test space of virtual displacements and the virtual interacting loads 
is the same space used for the approximations in Eqs. (49-66).) This process results 
in a 10M N + 13 M +1177+6 system of equations for the transverse shear deformation 
model and 6MN + 10 M + lOiV + 6 system of equations for the classical model, 
governing the displacements and the interacting loads. The governing equations 
are of the form 


Kn 

0 

0 

Bn 

B 12 

0 

A'22 

0 

B 21 

0 

0 

0 

K 33 

0 

B 32 

BTi 

bi ; 

0 

0 

0 

B ?2 

0 

Bl 2 

0 

0 

B \3 

Bl 3 

0 

0 

0 


B\3 " 


r a shell ' 


An >| 

B 23 


U$tr 


0 

0 

i 

Wring 

► = < 

0 

0 


^ str 


0 

0 


^ ring 


0 

0 . 


Q J 


. 0 , 


( 89 ) 
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in which sub-matrices A'n, A'22 and A' 33 are the stiffness matrices for the shell, 
stringer, and ring, respectively. The sub-matrices Bij,i,j = 1,2,3, in Eq. (89) are 
determined from the external virtual work terms involving the interacting loads, and 
the constraint Eqs. (46) to (48). The vector on the right-hand-side of Eq. (89) is the 
force vector, determined from the external virtual work terms involving pressure. 
The constraint equations correspond to the last three rows of the partitioned matrix 
in Eq. (89). Equation (89) is first solved for the displacements in terms of interacting 
loads, then this solution is substituted into the constraint equations to determine 
the interacting loads. Thus, the total solution is obtained. 

NUMERICAL EXAMPLE 

Numerical data for the example are R — 117.5 in., 21 = 20 in. and 2RQ = 5.8 in., 
which is typical of a large transport aircraft. The shell wall is a 16-ply quasi-isotropic 
[±45,0, 90, ±45,0, 90] s laminate of graphite/epoxy tape with a total thickness of 
0.080 in. The ply thickness is 0.005 in., and the lamina material properties are 
Ei = 1.85 x 10 7 lb/in. 2 , E 2 = 1.64 x 10 6 /6/m. 2 , G 12 = G 13 = 0.87 x 10 6 lb/in. 2 , 
G 23 = 0.49 x 10 6 lb/ in. 2 , and u i2 = 0.3. For the transverse shear deformation model, 
the shell wall stiffness sub-matrices of Eq. (23) are computed using these ply data 
and the expressions for the stiffness elements given in the Appendix. The numerical 
results are 

'0.654 0.198 0 

A = 0.198 0.654 0 x 10 6 Ibjin. 

0 0 0.228 _ 

'2.915 0 0.0904 0.0904' 

B = 0 -2.481 -0.0904 0.0904 lb 

0 -0.181 -0.652 0.652 _ 

■342.52 137.25 0 0 

137.25 291.53 0 -0.14 x 10~ 5 

D ~ 0 0 153.21 -0.30 x 10~ 5 

. 0 -0.14 x 10" 5 -0.30 x 10~ 5 0.30 x 10“ 5 

and Aij elements for the transverse shear stiffness matrix of Eq.(24) are 
A 44 — A55 = 0.544 x 10 5 Ibjin., A45 = 0 
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The bending and stretching-bending coupling sub-matrices for classical lamination 
theory, Eq. (35), are given by 


D 


342.52 137.25 

137.25 291.53 

0 0 


0 

0 

153.21 


lb in. 


B = 0 


The extensional stiffness sub-matrix A is the same for classical theory and shear 
deformation theory. Numerical data for the stiffeners are shown in Fig. 4. All the 
results presented are for an internal pressure p = 10 psi, and the Fourier Series are 
truncated at twenty- four terms in the x- and ^-directions (M = N — 24). 


RESULTS AND DISCUSSION 


Interacting Load Distributions 

The distributions of the interacting loads between the stringer and the shell 
are shown in Figs. 5 and 6. The distributions of the tangential component, A IS , 
shown in Fig. 5 are asymmetric about the origin (or stiffener intersection), and the 
maximum value in the transverse shear deformation theory is less by about 40% 
with respect to the maximum value in the classical theory. The distributions of 
the normal component, A zs , are also asymmetric about the origin (Fig. 6). The 
peak value of normal component computed from the classical model is 920.2 lb/in. 
and occurs at x = 0.2 inches. The peak value of normal component computed 
from the shear deformation model is 669.3 lb/in. and occurs at the origin. That 
is, the peak value of A Z3 is reduced by 27% in the transverse shear deformation 
model with respect to its peak value in the classical model. It is interesting to note 
that the asymmetric response predicted by the classical model is more predominant 
compared to the transverse shear deformation model. 

The distributions of the interacting loads between the ring and the shell are 
shown in Figs. 7 through 10. The distributions of axial force intensity, A rr , are 
symmetric about the origin, and attain extremum values at the origin as shown in 
Fig. 7. The extremum value of the axial force component is reduced from -144.3 
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lb/in. in the classical model to -20.3 lb/in. in the transverse shear deformation 
model (an 86% reduction). As shown in Fig. 8, the distributions of the tangential 
force component, A^ r , are antisymmetric about the origin, and the differences in the 
results for \e r from the two models are small. The distributions of the normal force 
intensity, X zr , are symmetric about the origin as shown in Fig. 9 and the predictions 
obtained from the two models are essentially the same. The normal force intensity 
is an extremum at the origin attaining a value of -2316.9 lb/in. in this example. 
The distributions of tangential moment component, Aq t , are symmetric about the 
origin, and attain extremum values at the origin as shown in Fig. 10. The maximum 
value of the tangential moment component is reduced from 525.2 lb-in. /in. in the 
classical model to 30.5 lb-in./in. in the transverse shear deformation model (a 
94% reduction). Results for \ xr and A 6r shown in Figs. 7 and 10, respectively, 
indicate that classical theory exaggerates the response of these variables, which are 
particularly sensitive to the asymmetry introduced by the ring. 

Stiffener Actions 

The distributions of the force and moment resultants in the stiffeners are shown 
in Figs. 11 through 13. The stringer axial force and bending moment distributions 
(Fig. 11) are slightly asymmetric about the origin, and only small differences are 
predicted between the two models. 

The distributions of the circumferential force and in-plane bending moment in 
the ring are shown in Fig. 12, and, again, there are only small differences in these 
results from the two models. However, the torque and out-of-plane bending moment 
in the ring are more sensitive to the change in models as shown Fig. 13. The 
distribution of the out-of-plane bending moment is symmetric about the origin and 
has reduced magnitudes in the transverse shear deformation model compared to the 
classical model. In the vicinity of origin, the transverse shear deformation model 
predicts less severe gradients in the distribution of out-of-plane bending moment 
compared to the classical model. As shown in Fig. 13, the distribution of torsion 
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is antisymmetric about the origin, and has reduced magnitudes in classical model 
compared to the shear deformation model. 

Shell Response 

The distribution of the normal displacement along x-curve midway between the 
stringers (0 = - 0), and along the 0-curve midway between the rings (x = - /), 
are shown in Fig. 14. As depicted in this figure, there is a negligible difference be- 
tween the results from the transverse shear deformation model and classical model. 
Also, there is negligible difference in the axial and circumferential normal strain 
distributions between the two models as shown in Figs. 15 and 16. 

A Ring with Symmetric Cross Section 

As a benchmark for comparing transverse shear deformation model with the 
classical model, analyses were performed for a ring with symmetric cross section. 
In this case the only change made to the example under discussion is to set the 
bending-coupling stiffeness EI ZX of the ring to zero from its value given in Fig. 
4. Consequently, the 0-axis, as well as the x-axis, are axes of symmetry for the 
repeating unit. Symmetry about the 0-axis implies there is no out-of-plane bending 
and torsion of the ring; i.e., u r (0) = <t>0r(@) — ^r(^) — ^ir(^) = A^ r (0) = 0 for 
-0 < 0 < 0 . 

The distributions of the tangential interacting loads between the shell and the 
stiffeners are shown in Figs. 17 and 18. For the tangential force intensity be- 
tween the stringer and shell (Fig. 17), the difference between the results from the 
transverse shear deformation model and classical model are small. For tangential 
force intensity between the ring and shell (Fig. 18), the differences in the results 
from the two models are larger than those for the stringer-shell tangential force 
intensity. However, differences between the two models for the ring-shell tangential 
force intensity are not excessively large. The tangential force intensity components 
between the stiffeners and shell are the most sensitive variables to a change in struc- 
tural model. All other dependent variables are only affected by negligible amounts 
to the change in structural models. 
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CONCLUSIONS 


The results for the shell-stringer interacting loads \ xs and \ zs , and shell-ring 
interacting loads \ xr and Ag r , indicate that the classical structural model exagger- 
ates the asymmetric response compared to the transverse shear deformation model. 
The predicted extremum values of these interacting load components are smaller in 
the transverse shear deformation theory than in the classical theory. The large dis- 
crepancy is due to the freedom in the transverse shear deformation model between 
the stringer’s bending rotation and the ring’s torsional rotation at the stiffener in- 
tersection, or joint. In the classical model these joint rotations are constrained to 
be the same. Shell strains and displacements are not significantly different in the 
two structural models. 
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Appendix 


ELEMENTS OF STIFFNESS MATRIX FOR A CYLINDRICAL SHELL 
BASED ON TRANSVERSE SHEAR DEFORMATION THEORY 

Based on the transverse shear deformation theory, the elements , and 

Dij of the stiffness matrices, in Eqs. (23) and (24) for the constitutive law for a 
laminated shell wall, are given by 
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#66 -^< 5662(1 + ^)(1 + ^) dz 
#66=^066^(1 + -^) <# 

#11 = J^Qnz 2 (l + 

0l2^ 2 ^ 

Qi6^ 2 (l + ^)dz 


#12 

£> 


D 2 

•^16 


#22 
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R' 

Qn~dz 

Q 2 1Z 2 {\ + dz 
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!e = /' 

-l 
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4^ (1 + S ) * 


where Qij are the transformed reduced stiffnesses given in the text by Jones (1975). 

Based on the assumption of constant transverse shear strain distribution through 
the thickness, the transverse shear stiffnesses are given by 

A 44 = Jc 44(1 + jj )^ 2 

J 45 — ^ C 45 J 2 

A55 = Jc5s(l + —) dz 


where 


C44 = G 1 3 C os 2 Q! -f" Gr2^Si'n?' OL 
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C 4 5 =(Gi 3 - G 2 3)CosaSma 

G55 =G 2 %Cos i 'OL + G 13 Sin 2 a 
in which a is the ply orientation angle. 

Based on the assumption of constant transverse shear stress distribution through 
the thickness, the transverse shear stiffnesses are given by 



in which k = knk 2 2 — k 2 2 - The coefficients kij are given by 

If z -1 

fcll J C 44(l + dz 

k\2 =J5 j C4$dz 

kn =£ /<*,( 1 + §)<fc 


where 


Cos 2 a Sin 2 a 
C44 — 7; H 


C45 =( 


(?13 

1 


1 


G13 G23 

Sin 2 a Cos 2 a 


G 23 

)CosaSina 


C55 


+ 


r 13 


<?23 
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Fig. 3. Displacements and rotations for (a) shell, (b) stringer, and (c) ring 
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Ring 

EA = 5.92 x 10 6 lb 
£/ xx = 26.9x 10 6 lb-in. 2 
EI^ = 2.69 x 10 6 lb-in. 2 
EI XZ = - 8.29 x 10 6 lb-in 2 
G7 = 1.0 x 10 6 lb-in. 2 
e r = 3.78 in. 


Fig. 4. Stiffener data. 


Stringer 

EA = 4.04 x 10 6 lb 
EIqq = 14.2 x 10 6 lb-in. 2 
(GA) X7 = 1.23 x 10 6 lb 
e s = 1.10 in. 



(Ib/in.) 





(Ib/in.) 





(Ib/in.) 





0/0 


Fig. 8 Ring-shell tangential force intensity. 
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Fig. 10 Ring-shell tangential moment intensity. 




Fig. 11 Stringer 
distribution. 
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Fig. 14 Distribution of shell's norm 
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Fig. 15 Distributions of the axial normal strain on the inner and 
outer shell surfaces midway between the stringers (0 = - 0). 




Fig. 16 Distributions of the circumferential normal strain on the 
inner and outer shell surfaces midway between the rings (x = - I). 
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Fig. 18 Ring-shell tangential force intensity for 
a ring with symmetrical cross section. 



